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The influence of intrinsic Fe defects in FeO (either single cation vacancies or prototypical 4:1 
vacancy clusters) on the electronic and dielectric properties is studied within the density functional 
theory. The importance of local Coulomb interactions at Fe atoms is highlighted and shown to 
be responsible for the observed insulating Mott gap in FeO which is reduced by the presence of 
defects. We investigate nonstoichiometric configurations of Fei-^O with x ranging from 3 to 9% 
and find the aliovalent Fe cations in both the regular and interstitial lattice sites of the considered 
configurations. Furthermore, we show that the trivalent Fe ions are induced by both isolated and 
clustered Fe-vacancies and introduce the empty band states inside the insulating gap, which de¬ 
creases monotonically with increasing cation vacancy concentration. The Fei-zO systems with high 
defect content become metallic for small values of the Coulomb interaction U, yielding the increase 
in the dielectric functions and optical reflectivity at low energies in agreement with the experimental 
data. Due to the crystal defects, the infrared-active transverse optic phonons split and distribute 
over a wide range of frequencies clarifying the origin of the exceptionally large spectral linewidths 
of the dielectric loss functions observed for wiistite in recent experiments. 

PACS numbers: 63.20.-e, 63.20.D-, 61.05.F-, 77.80.B- 


I. INTRODUCTION 

Wiistite, Fei-^O, has challenged experimental and 
theoretical condensed matter physics for over sixty years. 
In particular, the extent of its nonstoichiometry and its 
defect structure have been the subject of numerous exper¬ 
imental studies m and are still a matter of discussion 
and a source of controversy [6]. Generally, these stud¬ 
ies indicate that the crystal lattice of FeO is modified 
by the presence of Fe 3+ cations and vacant Fe 2+ sites. 
The crystal structure and charge distribution in Fei_ x O 
are very complex as the defects are likely to coalesce and 
arrange into clusters, which in turn can aggregate into 
larger defect structures. Although the size, shape and 
distribution of the stable equilibrium defect clusters and 
their aggregates are still controversial, they are expected 
to vary with pressure, temperature and composition [7], 
Moreover, the clusters of defects are believed to order 
magnetically below the Neel temperature Tjy = 198 K. 

Early theoretical investigations on nonstoichiometric 
FeO were mainly aimed at the mechanism of formation 
and stability of isolated Fe vacancies, prototypical 4:1- 
type clusters comprised of interstitial (tetrahedral) Fe 3+ 
ion surrounded by four cation vacancies as well as com¬ 
plexes of such clusters E9EU. The results of these studies 
indicate that wiistite favors formation of 4:1-type clusters 
which are major building blocks for further agglomera¬ 
tion and defect cluster growth [I]. Nevertheless, a sub¬ 
stantial concentration of free vacancies can also be stabi¬ 
lized within the wiistite lattice. 

Extended defect structures in Fei-^O were rarely mod¬ 


eled using the density functional theory (DFT). Such 
simulations are still extremely demanding, as they re¬ 
quire proper description of the highly correlated nature 
of 3 d electrons mm and extensive ionic relaxations 
within large supercells to predict a correct ground state 
of Fei_a;0. In our previous work, the DFT investigations 
have been performed for vacancy-defected FeO m , un¬ 
covering pronounced changes of the electronic and vibra¬ 
tional properties induced by isolated cation vacancies. 
Recently, it was shown that polaronic distributions of 
charge resembling those in magnetite Fea 04 emerge for 
the most stable defect structures in wiistite [IS] . 

The point defects as well as their clusters are likely to 
affect the dielectric and optical properties of FeO. The re¬ 
sults of various experiments [TTil - fITI indicate some pecu¬ 
liar behavior of the measured optical quantities in wiistite 
and assign it to the native defects. However, a detailed 
explanation of the influence of such defects on the elec¬ 
tronic and dielectric properties of FeO remains unclear. 
Furthermore, sensitivity of both the electronic structure 
and dielectric functions to the defect states as well as 
the extent of modifications induced by isolated or cluster 
defects still remains unexplored. 

The present work extends description of wiistite by 
investigating its electronic structure and dielectric prop¬ 
erties as functions of the concentration of cation mono¬ 
vacancies and vacancy clusters. Our theoretical studies 
take into account strong local interactions of Coulomb 
type to properly describe the electronic properties and to 
better understand the mechanisms of the defect-induced 
changes in Fei-^O. 
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The paper is organized as follows. The details of the 
calculation method are described in Sec. hd Section Mil 
characterizes electronic properties, whereas Sec. |IV A] and 
|IV B | present theoretical derivation of the dielectric prop¬ 
erties and comparison with the experimental data, re¬ 
spectively. Finally, Sec.[V]summarizes the current results 
and provides conclusions. 


II. METHODOLOGY 

Calculations were performed within the spin-polarized 
DFT method using the VASP code [25]. Electron-ion 
interactions were described in the framework of the 
projector-augmented wave (PAW) method [2B]. Valence 
electrons of Fe and O atoms were represented by the 
(3d 7 4s 1 ) and (2s 2 2p 4 ) configurations, respectively. The 
gradient-corrected exchange-correlation functional in the 
form proposed by Perdew and Wang (GGA-PW91) [27] 
together with a plane-wave expansion up to 520 eV 
were applied. Effects of electron correlations beyond the 
GGA were taken into account within the framework of 
GGA+t/ and the approach of Dudarev et al. [2S|, where 
a single parameter U e g = U — J determines an orbital- 
dependent correction to the DFT energy with J = 1 eV 
denoting the local exchange interaction. Most of the 
present calculations were performed using U e g = 5 eV, 
which is recognized as a realistic value while describing 
electron correlations between Fe(3d) states [121114] . Fur¬ 
thermore, it produces the lattice constant a = 4.35 A, 
which remains in a good agreement with various exper¬ 
imental studies [3®. Additional calculations with re¬ 
duced values of U e g , ranging from 2 to 4 eV, were car¬ 
ried out to investigate the influence of strong electron 
interactions and the screening processes due to defects 
on both the electronic structure and dielectric functions 
of wiistite. 

A variety of experiments performed over decades indi¬ 
cates that FeO is always deficient in iron and the con¬ 
centration of vacancies in the oxygen sublattice is sev¬ 
eral orders of magnitude smaller than the concentration 
of iron vacancies A. 2, Following the experimental ev¬ 
idence, we introduced vacancies into the cation sublat¬ 
tice of FeO and left its anion sublattice defect-free. The 
Fei_ x O structures with vacancy concentrations x = 3, 
5, 6, 9% were modeled by supercells [TH 125115U] , Each 
initial supercell of Fei-ajO was derived from the basic 64- 
atom FeO supercell which takes into consideration the 
antiferromagnetic (AFII) order JFlJ and comprises two 
ferromagnetic Fe-sublattices differing in the orientation 
of the spin magnetic moments (spin-up and spin-down 
Fe-sublattices). Neutral cation vacancies were created 
by removing the appropriate number of Fe atoms from 
the AFII supercell, whereas interstitial Fe atoms are in¬ 
serted into the empty tetrahedral positions of such super¬ 
cell. The AFII supercell lacking one Fe atom corresponds 
to Fei_ x O with x = 3%, while that with two Fe atoms 
absent conforms to Fei-^O with x = 6%. In the latter 
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FIG. 1. (Color online) Schematic representation of the proto¬ 
typical 4:1-type cluster with a single Fe ion in a tetrahedral 
configuration of Fei-^O. 


case, each of the ferromagnetic sublattices contains one 
Fe vacancy (Fpe)- There is, however, a number of pos¬ 
sible configurations within the simulated supercell with 
x = 6% differing among each other by VFe~ VFe distances. 
For further considerations we have selected the configu¬ 
ration minimizing the system energy and corresponding 
to the Vp e — Vp e distance of 1.22a. Details of such calcu¬ 
lations as well as more comprehensive discussion is given 
in our previous work m- 

The Fei-2,0 structures with x — 5% and x ~ 9% corre¬ 
spond to the 4:1-type clusters. The prototypical 4:1 clus¬ 
ter is schematically shown in Fig. [T] It is composed of an 
interstitial Fe atom surrounded by four cation vacancies 
arranged into a tetrahedron [5] . We note that due to the 
AFII magnetic ordering, three cation vacancies occupy 
the spin-up (or spin-down) ferromagnetic sublattice while 
one cation vacancy is always located in the ferromagnetic 
sublattice with the reversed spin direction. The intersti¬ 
tial Fe ion may take on the spin direction of either spin-up 
or spin-down cation sublattice. The AFII supercell con¬ 
taining one 4:1-type cluster corresponds to the overall 
cation deficiency of about 9%. Due to the cubic symme¬ 
try of the perfect FeO crystal, all possible locations of 
the 4:l-type cluster in the supercell are equivalent. The 
lower concentration of such defects, x = 5%, is simulated 
by placing one 4:1 cluster into the supercell elongated in 
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TABLE I. The list of considered Fei-zO configurations with 
the cation vacancy concentration (ir), types of defects, sizes 
of supercells with respect to the crystallographic unit cell of 
defect-free FeO, total number of atoms in the supercell (N), 
and number of Fe atoms (Np e ), 


X (%) 

type of defect 

supercell size 

N 

Np e 

0 

stoichiometric 

2x2x2 

64 

32 

3 

1 Upe 

2x2x2 

63 

31 

5 

4:1-type cluster 

4x2x2 

125 

61 

6 

2 Vp e 

2x2x2 

62 

30 

9 

4:1-type cluster 

2x2x2 

61 

29 


one of the main crystallographic axes. For convenience, 
we summarize the methodology of constructing isolated 
and clustered defects of wiistite in Table U 

Structural relaxations were performed for fixed vol¬ 
umes of the supercells followed from the optimized value 
of the defect-free FeO lattice. This approximation holds 
for the range of nonstoichiometries considered in the 
present work as the lattice constants of wiistite remain 
weakly dependent on x cam- The atomic positions were 
relaxed without symmetry constraints imposed. Only in 
the case where we determine the formation energies of 
defects, both volumes and atomic positions of the super¬ 
cells were relaxed. However, this procedure did not affect 
the final results as the lattice parameters of the super¬ 
cells containing defects changed by less than 1%. The 
Brillouin zones were sampled with the fc-point meshes of 
2 x 2 x 2 (x = 3, 6, 9,%) and 2 x 2 x 1 (x = 5%) gener¬ 
ated according to the Monkhorst- Pack scheme. Conver¬ 
gence criteria for the residual forces and total energies of 
particular systems were set to 0.01 eV/A and 0.1 meV, 
respectively. The complex dielectric functions as well as 
the intensities of the infrared-active phonon modes were 
obtained within the methodology proposed by Gajdos et 
al. [32] and implemented in the VASP code. 

The valence charges of Fe ions have been determined 
from the calculated electron contact densities at the Fe 
nuclei obtained within the full-potential (linearized) aug¬ 
mented plane-wave plus local orbitals [FP-(L)APW + lo] 
method implemented in the WIEN code [33] ■ The wave 
functions have been expanded into spherical harmonics 
inside the nonoverlapping atomic spheres having the radii 
of muffin-tin sphere Rmt and in the plane waves within 
the interstitial region. The Rmt belonged to the ranges 
of 1.93-2.03 a.u. for Fe and 1.66-1.74 a.u. for O atoms, 
with a value depending on the simulated system. The 
maximum l value for the expansion of the wave functions 
into the spherical harmonics inside the Rmt spheres was 
set to Imax = 10, while for the expansion of the wave func¬ 
tions within the interstitial region the plane-wave cutoff 
parameter K max = 7/ Rmt was applied. The charge den¬ 
sity was Fourier-expanded up to G max = 12 Ry 1 / 2 . Here, 
calculations were performed for the effective Coulomb 
interaction U e g = 5 eV, introduced in the rotationally 


invariant form proposed by Anisimov et al. 1341 . Calcu¬ 
lations have been carried out for the supercell volumes 
adopted from the pseudopotential method calculations. 
Only the atomic position were relaxed within the FP- 
LAPW methodology with the force convergence of 0.01 
mRy/a.u. and the symmetry constraints removed. 

The valence states of Fe ions in Fei-^O systems were 
identified according to the isomer shift (<5) systemat- 
ics [35] ■ In general, the isomer shift is expressed as 
6 = a(p — po), where p and po stand for the electron con¬ 
tact densities at the resonant nucleus in a given matrix 
and reference material, respectively. The symbol a de¬ 
notes the calibration constant characteristic for a partic¬ 
ular nuclear transition. Our calculations are performed 
for the 14.41-keV transition in 57 Fe and the isomer shifts 
are given with respect to metallic bcc a—Fe upon apply¬ 
ing previously determined a = 0.291 (a.u.) 3 mm/s [36]. 
More technical details of such calculations can be found 
in our earlier papers j30j |36j. 


III. ELECTRONIC STRUCTURE 

All considered Fei_ x O compositions reveal two differ¬ 
ent types of cation valence states, namely high-spin Fe 2+ 
and high-spin Fe 3+ residing both in the regular (R) and 
interstitial (I) positions of the cation sublattice. The 
compositions with 3 and 6% of isolated vacancies show 
both the Fe 2+ and Fe 3+ ions in R sites of the lattice, while 
those having 5 and 9% of vacancies and corresponding to 
the 4:1-type clusters exhibit Fe 2+ and Fe 3+ valence states 
in both the R and I lattice sites. Interstitial Fe becomes 
divalent when its spin is parallel to the spin direction of 
the magnetic sublattice from which three Fe atoms were 
removed, whereas it converts to the trivalent state if its 
spin is antiparallel to this sublattice (i.e., parallel to the 
ferromagnetic sublattice from which one Fe atom was re¬ 
moved). Thus, the interstitial Fe 2+ has the opposite spin 
direction to the interstitial Fe 3+ . 

The total number of trivalent ions results from the 
number of cation vacancies created in the lattice since 
each vacancy is accompanied by two trivalent config¬ 
urations. Therefore, six Fe 3+ ions occupy R sites of 
the wiistite cation sublattice containing 4:l-type clusters, 
with Fe 2+ ion at the I site and five Fe 3+ ions are located 
in the R lattice positions when the vacancy cluster in¬ 
volves interstitial Fe 3+ ion. The compositions x = 5% 
and x = 9% with interstitial Fe 3+ ions have lower en¬ 
ergies than the respective compositions with interstitial 
Fe 2+ ions by 14 and 8 meV/atom, respectively. This 
indicates that the wiistite lattice favors interstitial Fe 
cations in the trivalent charge state, while the intersti¬ 
tial Fe 2+ state can be considered as the metastable state. 
Nevertheless, it is interesting to note essential differences 
between the electronic structures of these two configura¬ 
tions that can be seen in Fig. [2] 

First of all, the unoccupied electronic states due to 
Fe 3+ ions appear in both the spin-up and spin-down 
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FIG. 2. (Color online) Densities of electron states calculated 
for Fei_xO containing 4:1-type vacancy clusters (x = 5%) at 
U e ff = 5 eV. Configuration with: (a) interstitial Fe 3+ ion, and 
(b) interstitial Fe 2+ ion. Positive (negative) values represent 
the spin-up (spin-down) components. The 3 d states of Fe 2+ 
and Fe 3+ ions residing in the regular (R) and interstitial (I) 
lattice sites are denoted by shaded areas (see the legend). 
Solid curve corresponds to 0 2 ~ states. The top of the valence 
band is taken as the reference energy (0 eV). 


channels for the system with the interstitial Fe 3+ ion 
[Fig- § a)]. On the other hand, these states exist only 
in one channel direction in the system with the intersti- 
tially located Fe 2+ ion [Fig. j2jb)]. In the latter case, the 
Fe 3+ ions induced in the R lattice positions occupy solely 
a single magnetic sublattice, whereas they are created in 
both magnetic sublattices in the former case. The in¬ 
terstitial Fe 2+ ions introduce their electronic states into 
the same energy bands as the Fe 2+ ions residing in the 
R lattice sites, as could be expected. Therefore, in the 
Fei-^O systems containing either 4:1 vacancy clusters or 
free vacancies [H] the bottom of the conduction band 
as well as the states lying at about —8 to —7 eV below 
the top of the valence band consist of states formed by 
trivalent Fe cations. The empty states in the conduction 
band located above ~ 3 eV from the top of the valence 
band arise from the divalent Fe cations occupying the 
R and I lattice sites. Each Fei-^O configuration shows 
significantly reduced insulating gap by the unoccupied 
electronic bands arising from the Fe 3+ ions stabilized by 
cation vacancies. 

Relative stability between monovacancies and vacancy 
clusters in Fei-^O has been determined from the respec¬ 


tive formation energies of neutral vacancies, defined as 
follows [37]. 


Eforin — Ed EygO T U Ep ei (1) 

where Ed is the total energy of Fei-^O supercell con¬ 
taining n vacancies, E Fe q is the total energy of the FeO 
supercell, and E Fe denotes the energy of Fe atom in the 
metallic Fe bcc crystal. It occurs that formation energies 
of systems with x = 3 and x = 6% are respectively higher 
by 107 and 36 meV per defect than the formation energy 
of system containing 4:1-type cluster (x = 9%). Thus, 
the wiistite lattice seems to prefer formation of vacancy 
clusters over the isolated cation vacancies. This result 
corresponds very well to the recent findings presented in 
Ref. [T5] . 

In principle, the divalent and trivalent Fe cations re¬ 
siding in the R and I sites of the wiistite lattice are 
distinguishable by methods sensitive to changes in the 
electronic structure of valence shell (e.g. the Mossbauer 
spectroscopy), which are a consequence of changes in 
charge/spin states, electron localization/delocalization, 
and the defect structure. Indeed, the results of our cal¬ 
culations given in Table [il] indicate stronger screening of 
s electrons by 3 d electrons for the high-spin Fe 2+ ions 
compared to the high-spin Fe 3+ ions, which results in re¬ 
spectively higher and lower values of the corresponding 
isomer shifts. Moreover, the isomer shift of the intersitial 
Fe 3+ is lower (0.384 mm/s) than the isomer shift of Fe 3+ 
occupying the R sites (0.487 mm/s). Similar behavior is 
encountered for Fe 2+ cations that exhibit slightly lower 
isomer shifts when located interstitially (1.039 mm/s) 
than octahedrally (1.108 mm/s). There is, however, no 
systematic dependence of 5(Fe 2+ ) on the vacancy con¬ 
centration x contrary to <5(Fe 3+ ) which slightly increases 
with increased x. 

The calculated isomer shift of the interstitial Fe 3+ ion 
in Fei_ x O is close to that measured for stoichiometric 
a-Fe 2 C >3 (0.37-0.38 mm/s) [55 as well as that of the 


TABLE II. Calculated (x = 0, 3, 6,9%) and experimental 
(bold numbers) [4] isomer shifts of the Fe 2+ and Fe 3+ cations 
occupying the R and I sites of wiistite. Values of the isomer 
shifts for cations in the R positions are averaged over respec¬ 
tive sites. Results are given for U e s = 5 eV. According to 
the 8 systematics [35], compounds containing high-spin Fe 2+ 
show <S in the range of 1.0-1.5 mm/s, while for the compounds 
with high-spin Fe 3+ the 8 ranges from 0.3 to 0.6 mm/s. 


X 

(%) 

<J(Fe 2+ ) 

(mm/s) 

5(Fe 3+ ) 

(mm/s) 

0 

1.057 (R) 

- 

2.9 

0 . 975 - 1.110 

0.45 

3 

1.103 (R) 

0.445 (R) 

4.8 

0 . 933 - 1.005 

0.60 

6 

1.097 (R) 

0.453 (R) 

9 

1.108 (R), 1.039 (I) 

0.476 (R), 0.382 (I) 
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TABLE III. The Bader charges and magnetic moments for 
different concentration of vacancies x. The values of magnetic 
moments obtained from the FP-LAPW method are given in 
parentheses. The Bader charges and magnetic moments are 
averaged over respective sites. 


monovacancies 


4:1 cluster 


X 

0 

3% 

6% 

5% 

9% 

Bader charges (e) 

o 2 - 

7.32 

7.28 

7.25 

7.28 

7.26 

Fe 2+ (R) 

6.68 

6.70 

6.71 

6.69 

6.69 

Fe 2+ (I) 

- 

- 

- 

6.62 

6.57 

Fe 3+ (R) 

- 

6.40 

6.38 

6.33 

6.33 

Fe 3+ (I) 

- 

- 

- 

6.30 

6.31 

Magnetic moments (/rs) 

Fe 2+ (R) 

3.74 (3.68) 

3.67 (3.60) 

3.69 (3.61) 

Fe 2+ (I) 

- 


- 

3.72 (3.64) 

Fe 3+ (R) 

- 

4.19 (4.17) 

4.20 (4.18) 

Fe 3+ (I) 

- 


- 

4.15 (4.12) 


tetrahedral Fe 3+ ion in stoichiometric Fe 3 C >4 at 4.2 K 
(0.39 mm/s) [35]. Furthermore, the isomer shift of Fe 3+ 
in the regular lattice sites of Fei_ x O is comparable to the 
experimental value of the octahedral Fe 3+ ion in Fe 3 C >4 
(0.51 mm/s). There is also a correspondence between 
(5(Fe 2+ ) in Fei_ x O and S of octahedral Fe 2+ in Fe 3 C> 4 , 
since the Mossbauer resonance spectra of Fe 3 C >4 obtained 
at 4.2 K reveal at least three components arising from the 
octahedral Fe 2+ ions, for which the isomer shifts range 
from 0.81 to 1.02 mm/s [35]. The sub-spectra due to the 
octahedral Fe 2+ ions results from the lowered symme¬ 
try of Fe 3 C >4 below the Verwey transition temperature 
Ty ss 120 K. We mention that the room temperature 
Mossbauer spectra of Fe 3 C >4 show the isomer shift for 
the octahedral cations of 0.66 mm/s resulting from the 
electron hopping between the octahedral Fe 2+ and Fe 3+ 
ions. It is advantageous to compare the calculated and 
measured isomer shifts of Fe cations in Fe-doped Coi_ x O 
which is also a compound from the family of the simple 
3d transition-metal monoxides, albeit showing much less 
degree of non-stoichiometry (1-3%) than wiistite. The 
high-spin Fe 2+ and Fe 3+ ions in Coi_ x O give rise to the 
Mdssbauer resonance lines having respectively the isomer 
shifts of 1.02 and 0.37 mm/s |30j. These values remain 
in close relationship with those measured by the emission 
Mossbauer spectroscopy (1.12-1.14 mm/s for high-spin 
Fe 2+ and 0.34-0.37 mm/s for high-spin Fe 3+ ) 40]. 

The Mossbauer spectra of realistic Fei_ x O samples 
01 U show much complexity as they are composed of 
several components reflecting quite a large number of 
different local environments accessible for Fe cations in 
the wiistite lattice. Although many attempts were un¬ 
dertaken [41], none of them described the local defect 
structure of wiistite in detail as the hyperfine parameters 
of the Fe cations deduced from the empirical fits could 


only provide a rough information on the statistically av¬ 
eraged environment of Fe 2+ and Fe 3+ in Fei_ x O. Thus, 
the Mossbauer spectra of samples with the same stoi¬ 
chiometry could be successfully fitted using quite distinct 
phenomenological models. A difficulty in the interpreta¬ 
tion of the experimental spectra presumably originates 
from a strong dependence of the wiistite defect structure 
on stoichiometry and the population of certain defect ar¬ 
rangements in the Fei_ x O lattice may vary with compo¬ 
sition in a complicated manner. Although, the present 
theoretical research considers a limited number of possi¬ 
ble defect configurations in wiistite, the results collected 
in Tabic [TT] can be useful while interpreting the Mossbauer 
spectra of Fei_ x O because they provide direct relation¬ 
ship between the local electronic structure of Fe cations 
(reflected by the respective isomer shifts), their location 
inside the wiistite lattice, and their immediate surround¬ 
ing. 


We have calculated the charges of ions using the Bader 
analysis (see Table III). Due to the hybridization effects, 
the obtained charges differ from the ionic valences and 
the difference between the Fe 2+ and Fe 3+ state is approx¬ 
imately 0.3 e. The trivalent cations exhibit enhanced spin 
magnetic moments comparing to those of divalent ones. 
We note that the magnetic moments remain dependent 
on U e ff but weakly depends on x. Within U e g ranging 
from 5 to 3 eV, they decrease linearly with decreasing 
U e g by 0.085 and 0.048 /is/eV for Fe 3+ and Fe 2+ cations, 
respectively. It is interesting to note that the magnetiza¬ 
tion isosurfaces around the Fe 2+ and Fe 3+ ions distinctly 
differ in their shapes, which can be seen in Fig. [3] The 
isosurface of Fe 3+ is almost spherically symmetric, while 
that of Fe 2+ contains cavities. The symmetric distribu¬ 
tion of magnetization reflects nearly equal occupation of 
five 3 d orbitals in the Fe 3+ state, whereas one additional 
electron in the orbital of Fe 2+ ion destroys such sym¬ 
metric population and contributes to the visible asym¬ 
metry in the magnetization isosurface. 


It is worth to analyze variation of the electronic prop¬ 
erties of wiistite as a function of the Coulomb interaction 
parameter U e g. Although the exact value of U e g remains 
unknown, it cannot be excluded that U e g may depend on 
the defect concentration due to the screening processes 
induced by delocalized charge carriers. The effect of U e g 
as well as the vacancy concentration on the band gap in 
Fei_ x O are displayed in Fig. [4] We observe that the en¬ 
ergy gap of each Fei_ x O composition decreases linearly 
upon decreasing U e g. Particular declines in the band 
gap energies are characterized by different slopes, i.e., the 
lower the x the higher the slope. Obviously, the struc¬ 
tures with higher vacancy concentrations exhibit more re¬ 
duced energy gaps comparing to that of defect-free FeO, 
mainly due to the increased number of empty Fe 3+ states 
induced in the gap. These states are located just below 
the conduction band dominated by the unoccupied elec¬ 
tronic states arising from the Fe 2+ ions that hybridize 
with the unoccupied states of oxygen anions (see Fig. [2]) . 
The best agreement with the experimental energy gap 
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(a) 


(b) 


FIG. 3. (Color online) Magnetization isosurfaces of Fe 2+ and 
Fe 2+ cations in Fei-^O containing 4:l-type vacancy clusters 
(a: = 9%) calculated at U e s = 5 eV. (a) Configuration with 
the Fe 2+ ion at the R site and the Fe 3+ ion at the I site, (b) 
configuration with the Fe 2+ ion at the I site and the Fe 3+ ion 
at the R site. The cation vacancies and the R sites of Fe 3+ ions 
are denoted respectively by white and red spheres, whereas 
magnetization isosurface is marked in green. For clarity the 
isosurfaces are drawn only for the I site and one selected R 
site. Note relatively random distribution of Fe 3+ ions around 
the 4:1 defect. 


of 1.15 eV, measured for Fei_ x O with x = 7% at 4.2 K 
[22] , is found for U e g = 4 eV. The unoccupied trivalent 
states constitute a band having the width increasing non- 
linearly with increasing vacancy concentration, as shown 
for typical U e g = 5 eV in the inset of Fig. [4j Simi¬ 
lar behavior is also observed for the remaining values of 
Ueg, although the lower values of the Hubbard interac¬ 



x(%) 


FIG. 4. (Color online) Influence of local interaction U e s and 
vacancy concentration x on the energy gap E g in Fei-zO. 
Inset shows the bandwidth (A) of the empty states of Fe 3+ 
cations as a function of increasing x for U e a = 5 eV. Dashed 
lines are guides for an eye. 


tion result in more pronounced broadening of that band. 
Application of U e g = 2 eV leads each Fei-^O composi¬ 
tion to a metallic state and to disappearance of trivalent 
states in the wiistite lattice. In this case, all Fe cations 
become divalent with the average value of the spin mag¬ 
netic moment of 3.57 /.is. 


IV. OPTICAL PROPERTIES 
A. Dielectric functions 

In general, the frequency-dependent dielectric function 
of a solid is a complex tensor, 

e(u) = e 1 (w)+i£ 2 (w). (2) 

It can be used to characterize the linear response of a sys¬ 
tem to an electromagnetic radiation. The real £i(w) and 
imaginary e 2 (w) parts of e(ui) describe dispersion and ab¬ 
sorption of the radiation in a given material, respectively. 
In addition, the zero-frequency limit of £i (ui) corresponds 
to the electronic part of the static dielectric constant of 
a material e<*,, a parameter being of fundamental impor¬ 
tance in many aspects of material properties, whereas 
£ 2 (w) is more specific and remains closely related to the 
band structure of the system. 

The calculated dielectric function Eq. ([2]) of wiistite 
may exhibit anisotropic behavior due to the rhombohe- 
dral distortion induced by the AFII ordering and the 
removal of symmetry constraints arising from the pres¬ 
ence of defects. Indeed, results of our calculations in¬ 
dicate subtle differences among the components of e(w) 
along the main crystallographic directions in the consid¬ 
ered vacancy-defected systems of FeO. Nevertheless, one 
can neglect this effect as too tiny to be observed exper- 
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imentally and consider the averaged values of £i(cj) and 
£2(w). 

Figure [5] shows that both components of e(u>) remain 
sensitive to the vacancy concentration in the regime of 
low photon energies (< 8 eV), while at higher energies 
the differences between dielectric functions in the systems 
with different x are insignificant. Although eoo increases 
with increased x, the shift upward does not exceed 10% 
within the studied range of vacancy concentrations. The 
calculated values of at U e g = 5 eV vary from 5.24 
to 5.53 — they correspond closely to the experimental 
£oo = 5.38 given by Hofmeister et al. m and theoret¬ 
ical value determined for perfectly stoichiometric FeO. 
On the other hand, one notes a broad range of val¬ 
ues deduced from a variety of experiments mum ns- 
123], which extends from 9.6 to 13. A considerably larger 
value of Coo in Fei-^O in comparison to typical ~ 5 
for the remaining simple transition metal oxides (MnO, 
CoO, NiO) used to be ascribed to a large Fe deficiency 
reaching up to several percent EH3 Much higher values 
of Coo derived from experimental spectra may indicate an 
existence of large amount of delocalized charge carries in 
the highly nonstoichiometric samples that can be easily 
polarized enhancing the value of £ 00 . One should also 
be aware that the static dielectric constant is measured 
at finite temperature and involves a contribution arising 
from phonons, whereas the majority of the calculations 
are actually performed for a static crystal (zero tempera¬ 
ture and neglecting zero-point vibrations) and hence they 
report values due to the purely electronic screening. This 
can be an additional source of errors, unless the phonon 
contribution is carefully extracted from the spectra. 

The £ 2 ( 0 ;) spectra clearly show the red shift of the 
absorption edge and reduction of the optical gap with 
increasing vacancy concentration. The decrease of the 
optical gap from ~ 2.2 to ~ 1.5 eV while going from 
defect-free FeO to Fei-ajO results from a presence of the 
Fe 3+ states in the band gap. A substantial reduction of 
the band gap in Fei_ x O with vacancy content ranging 
from 7 to 8% has been revealed by the recent infrared 
reflectivity and ellipsometry experiments [22] EH] • These 
studies provide value of the fundamental absorption edge 
in wiistite of 1.15 eV at 5 K and 1.0-1.3 eV at room 
temperature. 

The absorption spectra of Fei_ x O with different con¬ 
centrations of vacancies show many common features. 
First of all, the spectra are dominated by a wide and in¬ 
tense peak at 5.2 eV which shifts slightly to lower (higher) 
energies for Fei-^O systems containing isolated vacancies 
(vacancy clusters). Another similarity is connected with 
the broad and less intense peak at 6.7 eV moving upward 
by about 1 eV but only in wiistite containing vacancy 
clusters. Generally, both peaks diminish their intensities 
with the increased defect content. The absorption spec¬ 
tra exhibit broad distribution in the energy range of 10- 
25 eV and structureless plateau declining progressively 
at still higher photon energies. 

There are, however, prominent differences in £ 2 (u>) 
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FIG. 5. (Color online) Real and imaginary parts of dielectric 
function Q versus photon energy for (a) FeO with isolated 
vacancies, and (b) FeO with vacancy clusters. Calculations 
are performed with U e g = 5 eV. 


spectra between systems with isolated vacancies and dri- 
type clusters. They are encountered mainly below 4.5 eV 
and reflect differences in the electronic structures be¬ 
tween those systems (cf. Fig. [ 2 ] and Fig. 1 in Ref. [14] 1. 
In the absorption spectra of wiistite containing isolated 
vacancies the shoulder extending from about 1.5 to 4 eV 
originates from interband transitions that involve states 
of Fe 3+ and Fe 2+ ions residing in R lattice positions as 
well as the oxygen states that are always hybridized with 
cation states. The respective feature in wiistite with va¬ 
cancy clusters gets better resolved structure with a peak 
emerging at 4 eV and a tail with small swelling centered 
just above 2 eV. The 4 eV peak originates solely from 
the transitions due to the states of Fe 2+ ions at the R or 
I sites, provided the latter are present in the lattice [cf. 
densities of electronic states in Figs.[2](a) and |2j[b)] . 

On the other hand, the Fe 3+ states, either from R or 
I sites, that have the most considerable contribution in¬ 
side the band gap, give rise to transitions observed at 
energy in the vicinity of 2 eV. It is worth to note that 
the swelling at 2 eV immerses into the broad distribution 
when Fe 3+ ions occupy only the R sites. In such a case 
the low-energy £ 2 ( 0 ;) spectrum of wiistite with 4:1-type 
clusters resembles the £ 2 ( 0 ;) spectrum of wiistite with iso¬ 
lated vacancies. Therefore, the swelling can be assigned 
to transitions involving trivalent Fe cations at intersti¬ 
tial lattice positions. Moreover, an exact identification of 

























the transitions that are responsible for the peaks in the 
£ 2 (w) spectra of Fei-^O is hindered by a strong overlap 
and hybridization of the s, p, and d states in both the va¬ 
lence and conduction bands. This effect has already been 
pointed out for perfectly stoichiometric FeO by Rodl et 
al. [I3]. 


B. Comparison with experiments 

Optical absorption spectra extracted from the reflec¬ 
tivity measurements [2DJ on wustite with nonstoichiom¬ 
etry of the order of several percent show a very broad 
peak at about 4.5 eV and strongly enhanced absorp¬ 
tion at low photon energies (see Fig. [6]). The latter 
is absent in our calculated £ 2 ( 0 ;) spectra obtained for 
U e ff = 3—5 eV. Some discrepancy between our theoretical 
results for £2 (w) and those derived from the experiment 
m remains unclear especially that both £i(w) functions 
agree reasonably well in the line shape and magnitude 
within the energy range of 1-5 eV. We note, however, 
that the theoretical and experimental results are incon¬ 
sistent below ~ 1 eV due to abrupt growth of the experi¬ 
mental spectral amplitude of e(u>). Such peculiar behav¬ 
ior of both £i(u;) and £ 2 ( 0 ;) (not encountered in MnO, 
CoO, and NiO) was also observed in the data obtained 
from the ellipsometry spectroscopy [53] and has been at¬ 
tributed to defect absorption at Fe vacancies. Indeed, the 
experimental £i(w) and £2 (w) spectra exhibiting growth 
with oj —» 0 seem to reflect the narrow-gap nature of 
Fei_a;0 at vacancy concentrations of the order of several 
percent. 

The increase in absorption spactra at low energies, typ¬ 
ical for metals (Drude peak), indicates the existence of 
delocalized charge carriers despite of the Mott insulator 
state of wustite. These itinerant carriers (electrons or 
holes) exist due to the defected structure of Fei-^O. Such 
a complicated electronic structure is difficult to model 
within the DFT approach. However, the effect of wustite 
metallization can be to some extent simulated by lower¬ 
ing the Hubbard parameter U on the Fe cations. 

Results of such investigations are depicted in Fig. [6] 
One observes shifting of the peak’s positions to lower 
energies with decreasing U e s from 5 to 4 eV. At still 
lower U e g the spectra acquire more complex pattern com¬ 
posed of additional peaks appearing at low u>, and fi¬ 
nally the Fei_ x O system transforms to a metallic state 
at U e s = 2 eV . The £ 2 ( 0 ;) spectrum shows Drude-like 
peak at low energies in good agreement with the data 
derived from the experiment. The metallization arising 
from its large off-stoichiometry also affects the real part of 
the dielectric function which displays pronounced growth 
at low u, which in turn yields substantial increase of the 
electronic part of static dielectric constant (eoo=14.2). 
This phenomenon is likely to be responsible for the char¬ 
acteristic features of wustite dielectric functions detected 
in a variety of experiments and commonly interpreted as 
defect-related effects. 
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FIG. 6. (Color online) Dependence of (a) real and (b) imagi¬ 
nary parts of the dielectric function on U e s for Fei_a,0 with 
x = 6%. The values derived from the experiment (open sym¬ 
bols) are adopted from Ref. [20] . 


The real and imaginary components of e(ui) are usu¬ 
ally determined from ellipsometry experiments or de¬ 
rived from absorption, reflection, or transmission mea¬ 
surements after the Kramers-Kronig analysis. The ex¬ 
perimental data are frequently presented in the form of 
the reflectivity spectrum [42] . 


R(u) = 


- 1 ) 


RR + !) 


[n(uj) — 1\ + k 2 (uj) 
[n(w) + l] 2 + k 2 (oj)' 


(3) 


where 


»H = V , kMI + eiM . 
= y 


(4) 

(5) 


The functions ti(uj) and h(us) denote respectively refrac¬ 
tive index and extinction coefficients that constitute the 
complex refractive index, 


N{oj) = y/e(ui) = n(ui) + ik(ui). (6) 


A comparison between results of calculations per¬ 
formed as a function of x and the recent ellipsometry 
experiments carried out on samples with x = 7.5% [24] 
is given in Fig. [7] Similarly to the dielectric functions, 
the calculated spectra shift to lower energies for larger 



























9 



FIG. 7. (Color online) (a) Reflectivity R(uj ), (b) refractive 
index n(w) and extinction coefficient k(ui) of Fei-^O with dif¬ 
ferent nonstoichiometry. Theoretical and experimental data 
[24] are denoted by solid/dashed lines and symbols, respec¬ 
tively. Theoretical results are obtained for U e g = 5 eV. 

concentration of defects. Thus, the zero energy limit of 
the reflectivity and the refractive index increases and the 
absorption edge visible in the extinction coefficient di¬ 
minishes with the increasing x. Despite small differences, 
one finds that the calculated R{uj), n(u), and k(oj) spec¬ 
tra reasonably well reflect the behavior of the respective 
experimental quantities above ~ 1.5 eV. Both experi¬ 
mental and calculated data indicate an increase in the 
reflectivity of Fei-^O by almost 10% between 1.5-5.0 eV 
and the reflectivity decrease above this energy range. A 
visible mismatch between theoretical and experimental 
results is observed at low photon energies, where the mea¬ 
sured quantities show pronounced increase, probably due 
to higher amount of free charge carriers embedded in the 
measured samples as well as much stronger absorption 
by defects in comparison with the simulated Fei-^O sys¬ 
tems. 

A deficiency in the wustite cation sublattice was also 
suggested to be responsible for large spectral linewidth 
of the transverse optic (TO) phonon obtained from 
the infrared reflectivity measurements (22, 123]. One 
finds that the dielectric loss spectrum of the paramag¬ 
netic Fei_ x O contains broad and symmetric single line 
peaked at about 325 cm -1 , whereas the antiferromag¬ 
netic Fei_ x O reveals weakly resolved two-peak struc¬ 
ture making the low-temperature experimental spectrum 
asymmetric. This feature is also observed in our Gaus¬ 



200 220 240 260 280 300 320 340 360 380 400 420 440 
tOro (cm' 1 ) 

FIG. 8. (Color online) Theoretical infrared absorption intensi¬ 
ties (dotted vertical lines) of Fei-^O with the Fe-vacancy con¬ 
tent of 5% (vacancy clusters) and 6% (isolated vacancies) and 
the Gaussian convolutions of theoretical spectra with FWHM 
of 33 cm -1 (solid lines). Experimental dielectric loss function 
(221 measured at 5 K for Fei-ajO with x = 7 — 8% is repre¬ 
sented by solid symbols. The frequencies cjto correspond to 
the transverse optic phonons. Theoretical results are obtained 
for Ues = 5 eV. 


sian convoluted infrared absorption spectrum, which is 
directly related to dielectric loss function. Indeed, as 
depicted in Fig. [S] both theoretical and experimental in¬ 
frared absorption spectra show much broader distribu¬ 
tions of the TO phonon frequencies than those in more 
stoichiometric simple 3d transition metal oxides (MnO, 
CoO, and NiO) [23 . 

The broad theoretical spectrum arises from the split¬ 
ting of phonon frequencies due to the crystal symme¬ 
try being lowered by the presence of defects. The d- 
functions, shown in Fig. [8] represent the infrared in¬ 
tensity d-peaks associated with the splitted frequencies, 
subsequently convoluted with the experimental Gaussian 
width. It appears that the maximum of the convoluted 
spectral peak for the simulated compositions of x = 5 
and 6% shifts to higher frequencies (330 and 342 cm -1 ) 
in comparison with the position of the respective peak 
in the defect-free FeO (324 cm' 1 ). Also, the main peak 
in the Gaussian convoluted infrared absorption spectra 
changes its position from 338 cm -1 to 319 cm -1 , while 
increasing x from 3 to 9% (not shown). 

A similar trend to the one observed here has already 
been encountered in the early infrared reflectivity exper¬ 
iments m performed on wiistite samples with different 
concentrations of cation vacancies. Such a strong sensi¬ 
tivity of the infrared spectra to the concentration of de¬ 
fects can account for a large scattering seen in the exper¬ 
imentally determined TO phonon frequencies in wustite 
which cover the range extending from 320 to 410 cm' 1 
m IS [2TH23] . The broadening and shifting of the in¬ 
frared spectral peaks seem to have the same origin as 
very similar changes observed in the phonon density of 
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states M and they reflect pronounced modifications of 
the electronic structure and force constants due to cation 
vacancies in the Fe-sublattice of wiistite. 


V. SUMMARY AND CONCLUSIONS 

This work investigates the role of strong electron corre¬ 
lations on Fe atoms and the high concentration of cation 
vacancies in modifying the electronic and dielectric prop¬ 
erties of wiistite (Fei-^O). As we have shown, both of 
them influence substantially the electronic properties of 
this compound and have opposite effect on its band struc¬ 
ture. While the local electron interactions in the Fe 3 d 
states are responsible for the opening of the insulating 
gap, the Fe vacancies reduce its magnitude. The mecha¬ 
nism of gap reduction remains the same, irrespectively of 
the type of incorporated defects. Either monovacancies 
or vacancy clusters induce the band of empty Fe 3+ states 
inside the gap of the perfect FeO crystal. The width of 
this band increases monotonically with the increased con¬ 
centration of cation vacancies reducing the distance from 
the top of valence states, so effectively diminishing the 
band gap. 

The strength of the effective local interaction U e g de¬ 
pends presumably on the vacancy content, i.e., the band 
gap reduction enhances the screening of Coulomb interac¬ 
tions, which may lead to even larger electron mobility and 
metallization of iron oxide. The presence of free charge 
carriers, modeled here by small values of U e g, explain 
the increase of dielectric function and optical reflectivity 
observed experimentally at low energies. 

Summarizing, we presented an explanation of several 


effects observed in the optical properties of wiistite and 
explained them by the presence of defects. We empha¬ 
size that our studies also explain the anomalous broad¬ 
ening of infrared spectra, demonstrating a strong effect 
of cation vacancies on phonons in wiistite. Similar corre¬ 
lations between charge distribution and lattice dynamics 
induced by the electron-phonon coupling play an impor¬ 
tant role in the Verwey transition in magnetite (FeaCU) 
22]. As noticed recently also the polaronic distribu¬ 
tion found in defected FeO resembles the charge-orbital 
order in magnetite and may share a common origin with 
the short-range order observed above the Verwey tem¬ 
perature [22]. Finally, we point out that the effects dis¬ 
cussed here can be relevant not only for iron oxides but 
are expected to determine the electronic, dielectric, and 
dynamical properties of other strongly correlated systems 
with high concentrations of intrinsic defects as well. 
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